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Abstract 

The energy spectrum of eigenvalues of the equal spin-spin interac- 
tions (ESSI) Hamiltonian has been found. The obtained spectrum is 
free from limitations imposed on number of spins and parameters of 
the ESSI Hamiltonian. This model can be used for consideration of 
spin dynamics of mesoscopic systems and molecules with large number 
of nuclei spins. 

1 Introduction. Many-body problem of spin- 
spin interactions 

Thorough attention of researches to the properties of the Hamiltonian of 
spin-spin interactions is kept during some decades. It is due to the fact that 
this Hamiltonian determines fine peculiarities of dynamics and kinetics of the 
electron and nuclei spins of paramagnet spin-systems. The Hamiltonian of a 
paramagnet (electronic and nuclear) has the following form |]l| 

H = Hz + Hd + Held, (1) 

7iz = nuj^Y.^h '^^d = hJ2^fjS]S}, nM = hJ2^f,isrs+ + slsj), 
J fj fj 

= S''± iSy, Afj = Afj + Jfj, Afj = {h-i"^ I2R)-){1 - 3cos%), 
B;, = (-1/4)% + (1/2) 
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where Sj is the a-th spin component of a spin 5'^ (a = x, y, z), located in the 
point j, {Rfj, 9fj, ipfj) define the spherical coordinates of a vector connecting 
the points / and j in the laboratory frame which OZ axes is parallel to the 
applied constant magnetic field i^o, 7 is a gyromagnetic ratio, ojq = ■jHq is 
Zeeman frequency, the value Jfj defines an exchange integral. Only the sec- 
ular (diagonal) part of the dipole-dipole Hamiltonian is given by expression 
([^). This part determines the most essential contribution to dynamics and 
kinetics of spins in comparison with off-diagonal parts 0, |^, |^ . The Hamil- 
tonian (|ID represents a typical example of many-body interactions and so it 
is no needless to say about exact diagonalization of (|I]) if number of inter- 
acting spins n > 10. Even if the numerical calculation of the corresponding 
eigenvalues will be obtained, the analysis of this information becomes diffi- 
cult because the number of states = 2" is growing exponentially with the 
increasing of n. So in existing theories only macroscopic characteristics as 
spin-spin relaxation times, the second and the fourth moments of resonance 
lines were taken into account 0. 

Recently in investigation of the NMR on protons in a nematic liquid crys- 
tal a principal new NMR spectrum has been demonstrated. It consists of a 
set of well-resolved resonance lines instead of one usual unresolved resonance 
line broadened by spin-spin interactions. In recent time the authors brought 
the number of resonance lines up to the value 1024 0. This result opens new 
possibilities in studying of dynamics and kinetics of dipole-dipole connected 
spins in the mesoscopic systems and gives a new perspective physical medium 
for quantum information science. 

In such a way, the theoretical description of quantum-statistical properties 
of dipole-dipole coupled spin-systems becomes a very interesting problem. 

A liquid crystal where the spin-spin multiplet resolved spectrum with nar- 
row lines was observed |Q, includes partly oriented molecules, each molecule 
contains n = 19 interacting protons and intermolecular spin-spin interaction 
is negligible. It implies to find about a half-million of eigenvalues of the 
Hamiltonian that is the useless problem for exact solution. We want to 
stress that at certain values of orientation angles 9fj the dipole-dipole inter- 
action between close disposed spins can be less in comparison with remote 
spins. It is necessary to take also into account that orientation of spins in 
the space is not fixed precisely and true interaction constants are effectively 
replaced by their parts non-averaged by heat motion. This averaging of the 
spin-spin interaction of protons, located in the different part of the molecule 
is not equivalent: the effect of averaging is greater in the more mobile frag- 
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ments. So, the further progress is possible with using a simphfied theoretical 
model, which keeps nevertheless the main features of the basic Hamiltonian 
([1|). In our opinion such theoretical model can be presented by the equal 
spin-spin interactions (ESSI). The Hamiltonian of the ESSI has the form: 

n = i-Lz + i-Lad, [^) 

j fj fj 

where A and B defines the average values of parameters Afj and Bfj for 
separate molecule. The Hamiltonian of such type was considered in paper 
related to multiple nonergodic quantum dynamics . It seems that the range 
of applicability of this Hamiltonian is much wider viz., it can be applicable 
for consideration the sets of spins in mesoscopic magneto-active clusters. 

The averaged constants A and B take into account the interaction of the 
chosen spin with all others, keeping the difference between the longitudinal 
[A) and transverse [B) spin components. These constants can be obtained 
from real interaction potential as 

(A) = (1/n) 5^ A;„ AA = [{1/n) ^^(A;, - A^/^ A = (A) ± AA, 
fj fj 

and the same for constant B. 

In the next section the energy spectrum of the Hamiltonian (0) will be 
found. 

2 The matrix structure of the ESSI model 

The operators Tiz and Ti.^ commute with each other and with a longitu- 
dinal component Sj of a separate spin. So the eigenvalues 

Enp = nujo{2p - n)/2 + hA{3p'^ - 3np + - n) /4 (3) 

of the operator Hz + 'H'^ are expressed by means of eigenvalues (±1/2) of the 
operator S^, where n defines the total number of spins and p is a number of 
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spins oriented positively along the applied magnetic field. Eigenfunctions of 
the operator Ti^ + 

\Xa=i{n,p)) = 

= |mi = 1/2,1712 = 1/2, . . . , mp = 1/2, 
nip+i = -1/2, . . . , rrin-i = -1/2, m„ = -1/2) = 

= !+,+,•••,+,-, •••,-,-) (4) 

are also constructed as products of n eigenfunctions \mj) of operators Sj. 
The state (0) is degenerated z{n,p) = times at the given n and p. All set 
of eigenfunctions \Xa{n,p)), belonging to the given Enp is obtained from the 
function \Xa=i{n,p)) by permutation of its initial arguments running over all 
combinations from n over p, i.e. a = 1, 2, ... , z{n,p). 

Due to the structure of the operators + S'^Sj , the Hamiltonian 

does not have nonzero matrix elements between states with unequal number 
of positively oriented spins. So its matrix has blocked-diagonal form, where 
the sizes of blocks is determined as C^xC^. Diagonal elements of the operator 
Hz + 'H*^ inside blocks are the same and determined by expression (|^). In 
such a way the problem of diagonalization of the Hamiltonian (|^) is reduced 
to the finding of eigenvalues of the operator Ti*^^ in basis of states (H) with 
eigenvalues belonging to (Q). 

3 Eigenvalues of the operator TL*^^ 

For diagonalization of the operator Ti*^^ a special numerical program was 
written. It helps to find the desired matrices in basis (^) inside (?T,,p)-blocks 
and includes subsequent calculation of eigenvalues and eigenfunctions of cor- 
responding matrices. The spectrum of eigenvalues together with their de- 
generation multiplicity g of the operator has been found for all n and p 
satisfying to condition n < 8. Analyzing these results we established that in 
each (n,p)-block: 

1. the number of different eigenvalues is defined as 

r{n,p) = r{n,n — p) = p + 1, (5) 

2. the unequal eigenvalues of the operator are determined as 

hek{n,p) = hek{n,n — p) = hB[—p + k{n — 2p + 1) + A;^], (6) 
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3. the energy degeneration of the level ek{n,p) is 

9k{n,p) = gk{n,n-p) = C^-^ - Cl-^-\ (7) 

where < p < j)*, < A; < p, p* = [n/2\ is the integer part of n/2 and k 
numbers the operator Ti*^^ unequal eigenvalues. 

In further the correctness of these formulae was confirmed by direct nu- 
merical calculations for all possible values of n and p up to n = 11, p = 4. 
These calculations give us a ground to consider that formulae (j^), (|^) and 
(H) are exact expressions for the spectrum of the Hamiltonian describing 
ESSI-model. 

The eigenfunctions of the Hamiltonian 7i^^ were also calculated 

z 

\^p{n,p)) = ^C^p{n,p)\xa{n,p)), I < f3 < z{n,p), (8) 

a=l 

by determination the sets of coefficients 

{n, p) , C213 (n, p) , C^o (n, p) , . . . , C^/s {n,p)). 

Unfortunately, the compact expressions for coefficients Capiji^p) were not 
found. Expressions for Caf3{n,p) were obtained numerically for all cases, 
where the energy spectrum has been determined. These expressions are 
rather cumbersome and so, for example, in the Table 1 the coefficients 
Cai3in,p) were given only for n = 5. 

Taking into account the fact that the function Ixain^p)) belongs to the 
same eigenvalue (^, one can conclude that any their linear combination, 
in particular combination (|^), also belongs to this eigenvalue. Due to this 
property the eigenvalues of the ESSI Hamiltonian can be presented in the 
form of simple sum of energies (|^) and 

£k{n,p) = Enp + hek{n,p), (9) 

To each set of the three indices n,p, k corresponds gk{n,p) eigenfunction (H). 

In conclusion we would like to note that exact expressions (0), (|^)-(0) 
found for the quantum model ESSI are not exact solution but has the anzatz 
status. The found solution was guessed and then checked on large number 
of cases including block (n = ll,p = 4), which consists of 330 states. 
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Table 1: The eigenvalues, degeneration multiplicity value and eigenfunctions 
of the Hamiltonian for five spins 
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